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ABSTRACT 

Using the Chandra X-ray Observatory and the Rossi X-ray Timing Explorer, we have studied the black hole 
candidate (BHC) X-ray transient XTE J 1650-500 near the end of its 2001-2002 outburst after its transition to 
the low-hard state at X-ray luminosities down to L = 1.5 x lO-''* erg s"' (1-9 keV, assuming a source distance 
of 4 kpc). Our results include a characterization of the spectral and timing properties. At the lowest sampled 
luminosity, we used an 18 ks Chandra observation to measure the power spectrum at low frequencies. For the 
3 epochs at which we obtained ChandralRXTE observations, the 0.5-20 keV energy spectrum is consistent with 
a spectral model consisting of a power-law with interstellar absorption. We detect evolution in the power-law 
photon index from F = 1 .66 ± 0.05 to F = 1 .93 ± 0. 1 3 (90% confidence errors), indicating that the source softens 
at low luminosities. The power spectra are characterized by strong (20-35% fractional rms) band-limited noise, 
which we model as a zero-centered Lorentzian. Including results from an RXTE study of XTE J 1650-500 near 
the transition to the low-hard state by Kalemci et al. (2003), the half-width of the zero-centered Lorentzian 
(roughly where the band-limited noise cuts off) drops from 4 Hz at L = 7 x 10^^ erg s~' (1-9 keV, absorbed) to 
0.067 ±0.007 Hz at L = 9 X lO^^^ erg s"' to 0.0035 ±0.0010 Hz at the lowest luminosity. While the spectral and 
timing parameters evolve with luminosity, it is notable that the general shapes of the energy and power spectra 
remain the same, indicating that the source stays in the low-hard state. This implies that the X-ray emitting 
region of the system likely keeps the same overall structure, while the luminosity changes by a factor of 470. 
We discuss how these results may constrain theoretical black hole accretion models. 

Subject headings: accretion, accretion disks — black hole physics: general — stars: individual (XTE J1650- 
500) — stars: black holes — X-rays: stars 



1. INTRODUCTION 

To fully understand the behavior of black hole candidate 
(BHC) X-ray transients, it is important to determine the 
source properties for the full range of observed X-ray lumi- 
nosities from outburst peak at L ^ ]^q38-39 gj.g g-i quies- 
cence at L ~ 1 0^""^^ erg s"' . The brighter part of this range has 
been well-studied by several satellites, and, due to its broad 
bandpass (2-200 keV) and its excellent timing capabilities, 
the Rossi X-ray Timing Explorer (RXTE, Bradt, Rothschild 
& Swank 1993) has proved to be an extremely useful tool 
for studying the luminosity range above ^10"^^"^^ erg s"' (for 
sources with typical distances of 3-10 kpc). For most systems, 
this is below the luminosity where the source makes a tran- 
sition from the canonical high-soft state to the low-hard (or 
hard) state ( McClintock & Remillard 2003i) . and our group as 
well as others have used RXTE to study BHC transients as 
they undergo soft-to-hard state transitions during outburst de- 
cay (Tomsick & Kaaret 2000; Nowak. Wilms & Dove 200^ 
[Kalemci 2002; Kalemci et al. 2003a). However, as RXTE is 
not an imaging instrument, background and source confusion 
become significant at lower luminosities, making both spec- 
tral and timing studies problematic. Thus, at low luminosities, 
imaging observa tories such as Bepp oSAX, the Chandra X-ray 
Observatory ( Weisskopf et al. 2002), and XMM-Newton are 
necessary, and these sat ellites have been us ed to study BHC 
transients in quiescence JGarcia et alJl200ll) . The aim of this 
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Study is to use Chandra to study the evolution of XTE J 1650- 
500 in the luminosity range between the transition to the hard 
state and quiescence. 

Few previous studies have been dedicated to observing 
BHC transients at luminosities between the hard state tran- 
sition and quiescence even though the source evolution has 
not been well-characterized, and the physical changes that 
occur in this regime are unclear The results of the studies 
of BHC transients that have been undertaken in this regime 
are intriguing. A correlation between X-ray and radio flux 
has been fo und for several BHC systems in the hard state 
JCorbel et al. 2003: Gallo. Fender & Po olev 2003), and this 
has been interpreted as evidence for a connection between X - 
ray production and a compact radio jet fMarkoff et alJl2003l) . 
As the sources decay, it has been observed that the X-ray en- 
ergy spectra maintain a power-law shape; however, in some 
cases, there is evidence for spectral evolution with sources 
becoming softer at low l umin osities jE bisawa et al. 19941 
'Tomsi ck. Corbel & Kaarejl200lHCorbel et aL.2003^ .Kalemci 
2002). X-ray timing studies immediately after the hard state 
transition show that, in most cases, the characteristic frequen- 
cies drop ( Tomsick & Kaaret 2000 ; Nowak, Wilms & Dov9 
2002; Kalemci et al. 2003a|t), and this may indicate movement 
of the inner edge of the optically thick accretion disk away 
from the black hole. However, for the most part, the X-ray 
timing studies have not extended to luminosities much below 
the hard state transition. 

The hard state accretion geometry, the emission mech- 
anisms, and the connection between the accreting mate- 
rial and the compact jet outflow th at is presen t in the hard 
state for at least some BHCs (Fender 2001) are all cur- 
rently active areas of debate. Several theoretical mod- 
els have been suggested. One possibility is that a quasi- 



2 



J. A. Tomsick et al. 



spherical, optically thin region forms in the inner por- 
tion of the accretion flow as i n advection-dominated ac- 
cretio n flow (ADAF) models ("Nara van. McClintock & Yil 
The X-ray emission is predominantly due to ther- 
mal Comptonization in the hard state for ADAF models 
and for more generic "sphere-ndisk" models, and some suc- 
cess has b een achieved in comparing such models to ob- 
servations ( lEsin. McClintock & Naravan"1997^. 'Nowak et alJ 
[1999). However, for a range of disk viscosities, the ADAF 
model is convectively unstable ( O uataert & G ruzinov 2000), 
and this should be considered before concluding that any 
ADAF or sphere-ndisk model provides the correct physi- 
cal description. Another possibility that implies a much 
different site for X-ray production and a different accre- 
tion geometry is that the X-ray emission is due to mag- 
netic flares above the disk ( Galeev, Rosner & Vaiana 197^ 
|di Matteo, Celotti & Fabian 1999). Although the presence of 
a jet is not incorporated into the basic models described above, 
some recent wor k describes how a jet m ay arise from a "mag- 
netic corona" (Merloni & Fabian 2002). Other models have 
focused on the emission from the compact jet, and one sug- 
gestion is that the X-ray emis sion is produced in the jet vi a 
a synchrotron mechanism tMarkoff. Falcke & Fendeill2001l) . 
Improving the luminosity coverage for the hard state is im- 
portant for providing tests of these models. In addition, it is 
currently unclear whether the same physical model can apply 
to these systems in the hard state and in quiescence. Major 
physical changes to the system could be seen as another state 
transition at very low luminosities. On the other hand, the lack 
of a sharp transition could indicate that quiescence is simply 
a low luminosity extension of th e hard state as has been previ- 
ously suggested for GX 339-4 (ICorbel et alJ2000HKong et all 
I2?r00l) . 

We chose to carry out our program by observing XTE 
J 1650-500. This was a new BHC transient when it was 
discovered in outburst in 2001 (Remillard 2001). Although 
the compact object mass is not well-constrained so that it 
has not been confirmed that this is a black hole system, 
its X -ray properties strongly suggest th at it harbors a black 
hole (iTomsick et al."2003' Homan et al.! l2003t iKalemci et all 
|^03b). During its outburst, the source was optically iden- 
tified (Castro-Tirado et al. 2001), and a radio counterpart 
was found ( Groot et al. ,2001) . In quiescence, a radial ve- 
locity study was carried out, indicating a binary orbital 
period of 5.1 hr and an optical mass function of 0.64 ± 
0.03 M0 (Sanchez-Fernandez et al. 2002). The value for the 
mass function is relatively low for a black hole system and 
could indicate that the system contains a relatively low mass 
black hole, has a low binary inclination (/ < 40°), or both. 
We previously reported on RXTE observations of XTE J 1 650- 
500, including a study of its timing propert ies near the soft- 
to-hard state transition ( Kalemci et al. 2003b). We also re- 
ported on 14 day X-ray oscillations and ^100 s X-ray flares 
that were detected w ith RXTE near the end of the outburst 
llTomsick et al.l2003l) . The current work focuses on our stud- 
ies of the source with Chandra. 

2. OBSERVATIONS AND ANALYSIS 

We obtained 3 Chandra observations of XTE J 1650-500 
near the end of its 2001-2002 outburst, and the times of the 
observations are marked on the long-term X-ray light curve 
shown in Figure^ We compiled the 3-20 keV light curve in 
Figure [2 using all the observations of XTE J 1650-500 that 
were made with the Proportional Counter Array (PCA) on 




Fig. 1.— The 3-20 keV light cui-ve for the 2001-2002 XTE J1650-500 
outburst using PCA data from all the RXTE observations. The time of the 
soft-to-hard state transition, which occun'ed on MJD 52231.5, is marked with 
an arrow, and the gap in observations that is centered near MJD 52250 is due 
to solar angle constraints. After the gap, the source exhibited 14 day oscil- 
lations I Tomsick et al. 2003). The 3 Chandra observations occurred during 
these oscillations and are marked with squares. 



RXTE. The 3 Chandra observations were made 67, 78, and 
104 davs after the soft- to -hard state transition ( Homan et aO 
12003: Kalemci et al. 2003b), during the time when the source 
was exhibiting large-amplitude flux oscillations with a char- 
acteristic time scale of 14 days ( Tomsick et al. 2003). All 3 
observations were made at flux levels well below the level of 
the soft-to-hard state transition. 

The Chandra observations were made using the Advanced 
CCD Imaging Spectrometer (ACIS). In each case, the source 
was placed on one of the back-illuminated ACIS chips (S3). 
The first two observations (1 and 2), were made with the High 
Energy Transmission Grating (HETG) inserted, and each ob- 
servation lasted approximately 10 ks (see Tabled for exact 
times and other detailed information about the observations). 
To increase the throughput, the grati ng was not used for obser- 
vation 3, making photon pile-up ( Davisll20011) a concern. To 
mitigate the effects of pile-up, we reduced the CCD exposure 
time from 0.4 s per frame from the nominal 3.2 s per frame 
by using a 1/8 CCD sub-array. In addition, the telescope was 
pointed 2'. 7 away from the target, slightly blurring the point 
spread function (PSF) for XTE J 1650-500. An exposure time 
of 18 ks was obtained for observation 3. In addition to the 
Chandra observations, shorter 1.8-2.4 ks RXTE observations 
were co-ordinated to cover part of each Chandra observation. 

2.1. Chandra Analysis 

For the two grating observations, we used the "level 2" 
event list produced by the standard data processing (ASCDS 
v6.5.1 using CALDB v2.11). For both observations, we ap- 
plied a standard correction to the ACIS-S4 chip to correct a 
flaw in the chip readout. For this, as well as the data process- 
ing described below, we used the CIAO {Chandra Interactive 
Analysis of Observations) v2.3 software and CALDB (Cali- 
bration Data Base) v2.21. Prior to extracting the grating spec- 
tra, we produced background light curves. For observation 1 
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only, we found a few background flares where the count rate 
increased by 50-100% for a time period of less than 50 s. Due 
to their short duration, the flares make only a small contribu- 
tion to the total background. Thus, we did not remove these 
times in the subsequent analysis, and the full integration time 
was used for both observations. 

For observations 1 and 2, we produced 0.3-8 keV zero or- 
der images. We used the source detection algorithm wavde- 
tect to search for sources in a 2 by 2 arcminute region around 
the Chandra aim-point, and we only detect the target in each 
case. The observation 1 and 2 source positions are consis- 
tent with one another, and we obtain a mean source position 
of R.A. = 16h 50m OOs.98, Dec. = -A9° 57' 43".6 (equinox 
2000.0) with an uncertainty of 0".6. Th is position is con- 
sistent with the previously reported radio ( Groot et al. "2001!) 
and optical ( Castro-Tirado et al. 2001) positions, and we note 
that the Chandra position is somewhat more accurate than 
the 3" and 2" radio and optical positions, respectively. We 
conclude that the detected source is indeed XTE J 1650-500. 
Given the recent detections of X-ray jets from XTE J 1550- 
564 ( C orbel et al. 2002) and 4U 1755-338 (Anaelini & White 
l2003h . we examined the XTE J 1650-500 images for X-ray 
extensions in the vicinity of the source, but we do not see 
any features that are likely to be jets. Using the zero order 
ACIS response matrix and a spectral shape similar to that ob- 
served for the XTE J1550-564 jet (F = 1.6 with interstellar 
absorption), the upper limit on the 0.3-8 keV absorbed flux is 
6 X lO"''* erg cm'^ s"' at positions 10" from XTE J1650-500. 

We used the ClAO routine tgextract for the extraction of 
grating spectra and the ClAO routines mkgrmf and mkgarf 
to produce response matrices. Combining the +1 and -1 
MEG (Medium Energy Grating) 0.7-7 keV spectra yields 
4786 counts in 10001 seconds for observation 1 and 4342 
counts in 9509 seconds for observation 2. For the +1 and -1 
HEG (High Energy Grating) 1-9 keV spectra, there are 2706 
counts and 1872 counts for observations 1 and 2, respectively. 
We also produced estimated background spectra; however, we 
find that the background is not significant with estimates of 
-20 counts for the MEG and -1 1 counts for the HEG (0.4% 
of the detected counts in both cases). The higher MEG and 
HEG orders contain between —30 and —200 counts (source 
plus background). For each observation, the zeroth order con- 
tains —2500 counts, but there is a high level of photon pile- 
up. The flux obtained from the grating spectra implies that the 
number of counts would be a factor of 3 higher without pho- 
ton pile-up. Thus, due to the low number of counts in grating 
orders 2 and 3 and the high level of pile-up in the zeroth or- 
der, we focus on the order 1 spectra in this work. A final step 
in the extraction of the spectra is to apply a correction to the 
response matrix to account for the gradually changing low en- 
ergy ACIS response. We used the script corrarfXo apply this 
correction. 

For observation 3, we used the "level 2" event list produced 
by the standard data processing (ASCDS v6.6.0 and CALDB 
v2.12). After inspecting an image including the full energy 
band, we produced a 0.3-8 keV image to reduce the back- 
ground level. With the 1/8 sub-array and only the S3 chip 
active, the field of view is approximately 1' by 8'.4 (128 by 
1024 pixels), and we used wavdetect to search for sources. We 
found 4 sources detected at greater than A-a significance. The 
brightest source by far is at a position consistent with the XTE 
J1650-500 position found above for the grating observations. 
Although XTE J 1650-500 appears to be elongated, we used 
ClAO to simulate an off-axis PSF, and we found that XTE 



J1650-500 is consistent with a point-like source. The other 3 
sources also appear to be point-like, and the closest source is 
more than 1' from XTE J 1650-500. Despite using an instru- 
ment setup designed to mitigate pile-up, the measured ACIS 
count rate of 0.493 s"' for XTE J1650-500 (see Tabledl in- 
dicates that the data at the PSF core is moderately affected 
by photon pile-up. However, we were still able to use this 
data for spectral and timing analysis as described below and 
in Appendix A. 

2.2. Analysis 

We initially performed spectral and timing analysis for the 3 
RXTE observations that were made during the Chandra obser- 
vations. Using the "Standard 2" PCA data, we extracted en- 
ergy spectra using scripts developed at UC San Diego and the 
University of Tubingen that incorporate the standard software 
for RXTE data reduction (FTOOLS v5.2). We used the 2002 
February release of the "Faint Source" model for background 
subtraction. Obtaining the best possible background subtrac- 
tion is important for this work because of the low source count 
rates. To achieve this, we used Standard 2 data from the top 
anode layer only and excluded data from PCU 0, which has a 
higher background level due to the loss of its propane layer"*. 

In addition to the 3 RXTE observations, we extracted en- 
ergy spectra for the final 8 observations for which the PCA 
count rates are shown in Figure \l\ For these observations, 
the mean 5-8 keV PCA count rate is 0.021 ±0.018 s"' per 
Proportional Counter Unit (PCU) using all three anode lay- 
ers, indicating that the rate is consistent with zero at slightly 
more than the 1-cr level. We determined this rate in order 
to make a comparison to the rate expected from the Galac- 
tic ridge emission at t he Galactic latitude of XTE J 1650-500 
(b = -3.44°). Valinia & Mai'shall (1998) have determined the 
level of Galactic ridge emission (in PCA count rate) averaged 
over Galactic longitude as a function of b (see Figure 2 of 
Valinia & Marshall 1998). They found that the distribution 
could be described by narrow and broad gaussian compo- 
nents, and only the broad gaussian contributes at the Galac- 
tic latitude of XTE J 1650-500. Using the parameters for the 
broad gaussian leads to a best estimate of the 5-8 keV Galactic 
ridge rate of 0.71 s"' for all anode layers and 5 PCUs. Scal- 
ing the observed count rate given above (0.021 s"') to 5 PCUs 
gives 0.11 ±0.09 s"', which is significantly below the best 
estimate for the Galactic ridge emission. Thus, we conclude 
that the Galactic ridge emission can easily explain the low 
level of emission we detect for the final 8 observations. As 
this analysis shows that it is likely that XTE J 1650-500 does 
not contribute to the flux detected during these observations, 
we produced an average spectrum (top anode layer only) for 
the final 8 observations and subtracted it from the PCA spectra 
for the 3 RXTE observations that were simultaneous with the 
Chandra observations. The effect to the continuum is small; 
however, the Galactic ridge emission contains a strong iron 
emission line at 6.7 keV, and this line is apparent in the aver- 
age spectrum for the final 8 observations. 

Finally, we analyzed the data from the 15 RXTE observa- 
tions taken between MJD 52265 and MJD 52298, which is 
after the gap in RXTE coverage and before the Chandra ob- 
servations began (see FigureQ. We processed the PCA data 
as described above, including the subtraction of the Galactic 
ridge emission from the energy spectrum. For the spectral 

See http:/Aheawww.gsfc.nasa.gov/users/craigm/pca-bkg/bkg-users. 
for a detailed analysis of the background model performance. 



4 



J. A. Tomsick et al. 





Mi 

. . 1 




11 1 1 1 

_ ^JLip: 


1 1 1 1 1 1 1 1 1 : 





0.5 1 2 5 10 20 

Energy (keV) 



Fig. 2. — The observation 1 Chandra and PCA energy spectra fitted 
with an absorbed power-law model. The Chandra instruments used are the 
ACIS/HETG/MEG (0.7-7 keV) and the ACIS/HETG/HEG (1-9 keV). The 
measured power-law photon index is F = 1.66 ±0.05. The bottom panel 
shows the residuals. 
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Fig. 3. — The observation 3 Chandra and PCA energy spectra fitted with 
an absorbed power-law model. We obtained the Chandra spectrum (0.5-8 
keV) with direct ACIS imaging. The measured power-law photon index is 
F = 1 .93 it 0. 13. The bottom panel shows the residuals. 



analysis described in § 3.3, we also removed flares fr om 2 of 
the observations as described in lTomsick et alJ (120031) . 

3. RESULTS 
3.1. Energy Spectra at Low Flux Levels 

Using XSPEC 11.2, we performed simultaneous least- 
squares fits to the Chandra and PCA spectra for all three ob- 
servations. For the PCA, we included 0.6% and 0.3% sys- 
tematic errors for energy bins below and above 8 keV, re- 
spectively, as described in our previous work ( Tomsick et al. 
120031) . We began by fitting the observation 1 spectrum with 
a power-law with interstellar absorption. We also included 
a constant factor to allow for the possibility that the overall 
normalization is different for the PCA and Chandra, but we 
assumed that the overall MEG and HEG normalizations agree 
with each other As shown in Tabled a power-law model 
with r = 1.66 ±0.05 and A^h = (6.7 ±0.5) x lO^' cm'^ (er- 
rors are 90% confidence, Ax^ = 2.7) provides a good fit to 
the spectrum (x^/ v= 116/ 155). The other continuum models 
we tried do not provide acceptable fits. Thermal model s, such 
as a blackbody or a disk-blackbody (^Makishima et al.<ll986h . 
result in reduced values of 8.4 and 2.3, respectively. Fig- 
urelJlshows the observation 1 spectrum and residuals with the 
power-law fit. No clear features are seen in the residuals, and 
we note that the fit is not improved by adding a soft compo- 
nent, a high energy cutoff, or an iron emission line. The upper 
limit on the equivalent width of a narrow (where "narrow" is 
set by the binning shown in Figure |2} Fe Ka emission line at 
6.4 keV is 161 eV. 

Tableland Figure ^ indicate that the count rates are very 
similar for observations 1 and 2. We fitted the observation 2 
spectrum, and Table |2] shows that the spectral parameters do 
not change significantly between the two observations. Also, 
in both cases, the spectral fits indicate that the PCA/Chandra 
normalization factor is near 1.12. As the PCA exposure only 
covers part of the Chandra observation, we checked to see if 
this difference could be due to source variability. We extracted 
a grating spectrum for observation 1 that is strictly simultane- 
ous with the PCA exposure. During this time, the MEG count 
rate is 0.483 ±0.016 s"' (0.7-7 keV), which is consistent with 
the rate for the full observation (0.479 ±0.007 s"'). Thus, 
we conclude that the overall normalizations are different for 



the PCA and the Chandra gratings, and we note that this is 
not surprising since it has been previously determined that the 
PCA norma lization is high relative to previous instruments 
dTomsicket al..l999,) . 

Due to the presence of photon pile-up in the PSF core, 
we extracted the observation 3 ACIS spectrum after remov- 
ing the counts detected in the PSF core. We extracted the 
source counts from an elliptical annulus centered on the XTE 
J 1650-500 position. We used an inner ellipse with a semi- 
major axis of 1".3 and a semi-minor axis of 0".9, and an outer 
ellipse with a semi-major axis of 4". 9 and a semi-minor axis 
of 3".l. We rotated the ellipses to match the off-axis Chan- 
dra PSF. This region contains 1573 counts compared to the 
9000 counts obtained when the core is included (see Tabled, 
and we are confident that the photons from the elliptical an- 
nulus do not suffer from pile-up. However, it is critical to 
realize that ACIS response matrices produced by the standard 
ClAO software do not account for the size or shape of the 
source extraction region or the fact that the Chandra PSF is 
energy dependent. We corrected the ACIS response matrix for 
this using the Chandra Ray Tracing software (ChaRT) as de- 
scribed in Appendix B. When producing the response matrix 
we also applied the low energy ACIS response correction for 
detector evolution as described above for the grating spectra. 

We fitted the observation 3 ACISh-PCA spectrum with an 
absorbed power-law model as shown in Figure^ This model 
provides a good fit to the data (x^/ // = 104/100) with parame- 
ter values ofA^H = (6.4±0.7)x 10^' cm-2 and T = 1.93 ±0.13 
(see TableO. The measured column density is consistent with 
the value obtained from the grating observations. Thus, we 
re-fitted the spectrum after fixing the column density to the 
value found using the grating spectra (A^h = 6.7 x 10^' cm"^), 
and we obtain F = 1.96 ±0.09. Also, for consistency with 
observations 1 and 2, we used the PCA normalization for ob- 
servation 3 divided by the PCA/grating normalization factor 
of 1 . 12 to determine the observation 3 flux reported in Table|5] 
The spectral results indicate that the XTE J 1650-500 keeps its 
power-law shape but softens somewhat at low luminosity. 

3.2. Power Spectra at Low Flux Levels 

We examined the PCA and Chandra light curves for all 
three observations. As shown in Figure|3 clear variability is 
present in the PCA light curves for observations 1 and 2 and 
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Fig. 4. — PCA (panels a and b) and ACIS (panel c) light cui-ves for the 3 
observations. The energy bands are 3-20 keV and 0.5-8 keV for the PCA and 
ACIS, respectively, and the light curves have been background subtracted. In 
each case, a representative en'or bar is shown. 
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Fig . 5 . — (a) PCA power spectrum for observations 1 and 2 combined. The 
solid line shows the best fit zero-centered Lorentzian, and the power spectrum 
shows marginal (between 2 and 3-cr) evidence for a QPO at 0. 1 1 Hz. (b) ACIS 
power spectrum for observation 3 with the best fit zero-centered Lorentzian. 
In panel a, the dashed line shows the model used for observation 3 to illustrate 
that the power shifts to lower frequency for observation 3. 



in the ACIS light curve for observation 3. The plotted 3-20 
keV PCA light curves (panels a and b) have 4 s time bins, and 
we note that the count rate is higher than reported in Tabled 
because we used all anode layers and PCUs 0, 2, and 3 rather 
than just the top layer and PCUs 2 and 3. While the quality 
of the background estimate is slightly worse when PCU and 
all layers are included, the higher count rate allows for tighter 
constraints on the timing properties. We also produced grating 
light curves (MEGh-HEG) for observations 1 and 2. Although 
the statistical quality of the PCA light curves is superior to 
that of the grating light curves, we use the grating light curves 
for the analysis below as they allow us to study the energy de- 
pendence of the timing properties. Figure]^ shows the 0.5-8 
keV ACIS light curve for observation 3 with ~50 s time bins. 
We also examined the PCA light curve for observation 3, but 
it is background dominated and is not useful for timing anal- 
ysis. 

For observations 1 and 2, we produced power spectra us- 
ing PCA data taken in a high time resolution "event" mode 
(labeled E_125us_64M_0_ls). These data consist of an event 
list with photons time tagged to within 125 /xs and 64 en- 
ergy bins. First, we used the data in the 3-20 keV band to 
produce Leahy normalized power spectra (iLeahv et al.ll983|) 
for observations 1 and 2 separately. In each case, we aver- 
aged power spectra from four 460 s light curve segments and 
used a Nyquist frequency of 64 Hz. No noise in excess of 
the Poisson level is present above ~1 Hz, and there were 
no clearly significant differences between the power spec- 
tra for the two observations. Thus, we produced another 
Leahy normaUzed power spectrum (again with a 3-20 keV 



energy band) consisting of the eight 460 s light curve seg- 
ments from observations 1 and 2 with a Nyquist frequency of 
8 Hz. We converted to the rms normalization ( Mivamoto et aQ 
119941) . and this power spectrum is shown in Figure The 
power spectrum is dominated by a band-limited noise compo- 
nent, and such a component has typically been fitted with a 
broken power-law or a zero-centered Lorentzian in previous 
work (Belloni. Psaltis & van der Klis 2002). A zero-centered 
Lorentzian with a half- width (FWHM/2) of 0.067 ±0.007 Hz, 
and a fractional rms amplitude of 33.1%±1.1% (see Ta- 
ble |3} provides a reasonably good fit to the power spectrum 
(x^/f^ = 93/72). We note that the half-width of a zero- 
centered Lorentzian is equivalent to the "peak" frequency 
that is s ometimes quoted when fitti ng Lorent zians to power 
spectra JBelloni. Psaltis & van der Klis. .2002t iKalemci et alJ 
2003b). Positive residuals are present near 0.1 Hz, possibly 
indicating the presence of a quasi-periodic oscillation (QPO). 
The addition of a second Lorentzian at 0. 1 1 Hz improves the 
quality of the fit, resulting in x^/j^ = 83 /69; however, this im- 
provement corresponds to a detection at only the 96% confi- 
dence level. Although we repeated this analysis several times 
with different frequency binnings, we consistently found that 
the QPO is only detected at a level of between 2 and 3-(j. 
Thus, we take the measured rms amplitude of 12% as an up- 
per limit on a possible QPO at 0. 1 1 Hz. 

Using the grating light curves, we produced a 0.7-9 keV 
power spectrum for observations 1 and 2 by averaging eight 
2414 s light curve segments. The time resolution is set by 
the interval of time between chip read-outs, and this inter- 
val is 2.54104 s for these observations, giving a Nyquist fre- 
quency of 0.2 Hz. Although the time resolution and the sta- 
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tistical quality are lower for the grating light curves than for 
the PCA data, using the grating data allows us to extend the 
timing study to lower energies and frequencies. After convert- 
ing to rms normalization, we fitted the grating power spec- 
trum with a zero-centered Lorentzian, and the fit is accept- 
able ix^/v = 25/23). As shown in Table E] the half-width 
is 0.079 ± 0.020 Hz, and the fractional rms is 30.1%±2.5%. 
Thus, the parameters for the 0.7-9 keV power spectrum are 
consistent with those for the 3-20 keV PCA power spectrum. 

We also produced a power spectrum using the observation 
3 ACIS data. Using a 0.5-8 keV light curve with a time reso- 
lution of 0.44104 s, we produced a Leahy normalized power 
spectrum with a Nyquist frequency of 1.13 Hz, consisting of 
the average of six 3350 s light curve intervals. In the Leahy 
normalization, Poisson noise leads to a power level of 2.0; 
however, at high frequencies (0. 1-1 Hz) the ACIS power spec- 
trum is consistent with a constant level of 1 .77. We performed 
simulations described in Appendix A showing that pile-up at 
the level present for observation 3 leads to a drop in the Pois- 
son noise level from 2.0 to 1.78. This, along with a small 
drop in the Poisson noise level caused by deadtime during 
chip read-outs (see Appendix A) explains our measurement of 
1.77. Thus, when converting the ACIS power spectrum from 
the Leahy normalization to the rms normalization, we sub- 
tracted a value of 1.77 rather than the usual value of 2.0. The 
rms normalized power spectrum fitted with a zero-centered 
Lorentzian is shown in Figure [Sj", and this model provides 
a good fit to the data (y^ jv = 56/65). The Lorentzian half- 
width is 0.0035 ±0.0010 Hz (see Tableg), which is signifi- 
cantly lower than the value obtained for observations 1 and 2. 
This indicates that, as the flux drops, the band-limited noise 
moves to lower frequencies. This fact is illustrated in Fig- 
ure |5^, where the model used to fit the observation 3 power 
spectrum is plotted as a dashed line on the power spectrum for 
observations 1 and 2. The fit to the observation 3 power spec- 
trum indicates a fractional rms amplitude of 17.7%±1.7%, 
but this (16%) should be taken as a lower limit on the ac- 
tual value since the noise level is reduced by pile-up. From 
our simulations, we estimate that the actual rms amplitude is 
near 19%. We also produced a power spectrum using only the 
1573 counts from the region of the CCD that is free of pile-up. 
Fitting the power spectrum with a zero-centered Lorentzian, 
the half-width is not well-constrained, but it is consistent with 
the previous fit, and the rms amplitude is 19%±7%, giving an 
upper limit of 26%. Thus, considering the fits to the power 
spectra with and without pile-up, we conclude that the rms is 
19%!™. 

3.3. Summary ofXTE Jl 650-500 Hard State X-Ray 
Observations 

For observations 1-3, two of our main results are that the 
source softens (F changes from 1.66±0.05 to 1.93 ±0.13) 
and the band-limited noise component moves to lower fre- 
quencies (the Lorentzian half-width, FWHM/2, changes from 
0.067 ±0.007 Hz to 0.0035 ± 0.0010 Hz) as the source flux 
decays. Thus, where possible, we determined the values of 
these parameters (F and FWHM/2) for other RXTE obser- 
vations of XTE J 1650-500 in the hard state during outburst 
decay. For the 7 observations made immediately after the 
soft-to-hard state transition on MJD 52231.5 (see Figure 
we use the spectral parameters from Kalemci (2002). These 
energy spectra are dominated by a power-law component, 
but it should be noted that a soft component is also present. 
For these 7 observations, the timing parameters come from 



iKalemci etanj2003bh . In contrast to our ChandralRXTE ob- 
servations at low flux levels, several Lorentzians are required 
to describe the power spectra for most of the observations im- 
mediately after the soft-to-hard state transition. The lowest 
frequency component in those power spectra are nearly zero- 
centered Lorentzians, and we assume that they correspond to 
the band-limited noise component that we see in the low flux 
power spectra. 

For the 15 RXTE observations after the gap in RXTE cov- 
erage and before the Chandra observations began, the 3-20 
keV energy spectra are all well-described by a single power- 
law component with the column density fixed to the value of 
A^H = 6.7 X 10^' cm"^ found for the grating observations. For 
these observations, the 3-20 keV flux varies from 1.2 x 10~" 
to 2.9 X 10"'" erg cm"^ s"', and the values of F range from 
1.69 ±0.02 to ^2.0. The data suggest a trend of softening 
as the flux drops consistent with that observed with Chandra, 
but the errors on F are large for the low flux observations. 
We produced power spectra for all 15 observations; however, 
we could not obtain reliable fit parameters for about half of 
the observations due to poor statistics, short observations, or 
flaring (Tomsick et al. 2003). In 8 cases, we obtained use- 
ful power spectra, and these are dominated by a band-limited 
component with a fractional rms amplitude near 30%. We 
fitted these power spectra with a zero-centered Lorentzian to 
determine the half-width for comparison to the other observa- 
tions. Figure|6^ shows the measured values of the Lorentzian 
half- width vs. 3-9 keV flux, and Figure |6j3 shows the half- 
width vs. the measured values of F. 

For the data points above a half-width of 0.5 Hz, there is 
a tight correlation between F and the half-width with a linear 
correlation coefficient of 0.96. The two points near 0.25 Hz 
also appear to be consistent with this correlation, but the cor- 
relation coefficient drops slightly to 0.91 when these points 
are included. It is clear that the correlation does not extend to 
lower frequencies as F reaches a minimum near 1.65-1.70. In 
fact, the data at frequencies between 0.03 Hz and 0.3 Hz sug- 
gest the possibility of an anti-correlation between F and the 
half-width, but the linear correlation coefficient is only -0.41. 
We checked to see if it is valid to compare the ChandralRXTE 
point in this frequency range to the RXTE only points by fit- 
ting the RXTE spectra alone for observations 1 and 2. Fixing 
A^H to 6.7 X 10^' cm"^, we obtained values of F consistent with 
those found using Chandra and RXTE together Although 
the anti-correlation is not statistically significant, the distri- 
bution of the values of F for the 0.03-0.3 Hz points is also 
not well-described by a constant at their weighted average of 
F = 1.70 ±0.01. Thus, we conclude that there is some scatter 
in F for the 0.03-0.3 Hz points, but we cannot establish a def- 
inite relationship between F and the Lorentzian half-width in 
this frequency regime. 

4. DISCUSSION 

One of the main motivations for this work is to measure 
the X-ray properties of an accreting black hole in the rela- 
tively unexplored luminosity range between the transition to 
the low-hard (or hard) state and the quiescent state. Such mea- 
surements allow us to address the question of whether these 
should be considered as two separate states or if the quiescent 
state is, in fact, a low luminosity version of the hard state. This 
question has important implications for whether there are sig- 
nificant physical changes in accreting black hole systems that 
occur at low luminosities such as changes in the geometry of 
the accretion flow and/or emission mechanisms. 
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Fig. 6. — Measured pai'ameters from Chandra and RXTE observations of 
XTE J1650-500 in the hard state, (a) The values of the Lorentzian half-width 
vs. the flux in the 3-9 keV energy band, (b) The half-width vs. the power-law 
photon index (F). In both panels, the diamonds mark the measurements froin 
observations with only RXTE data and the filled circles mark observations 
where we obtained both Chandra and RXTE data (one point for observations 
1-1-2 and one point for observation 3). The eiTor bars on frequency and on F 
are 90% confidence. 



For XTE J 1650-500, we have measured the spectral and 
timing properties of the source at 1-9 keV (absorbed) lu- 
minosities of 9 X lO^'* erg s~' (observations 1 and 2) and 
1 .5 X lO-'"^ erg s"' (observation 3). The quoted luminosities as- 
sume a fiducial source distance of 4 kpc, and we note that the 
u nabsorbed luminositie s are about 25% higher As discussed 
in lTomsick eTaP (l2003h . the likely distance for XTE J 1650- 
500 is between 2 kpc and 6 kpc. For the RXTE observation 
immediately after the soft-to-hard state transition, we used the 
spectral model described in Kalemci (2002) to extrapolate the 
PCA-measured flux into the 1-9 keV energy band, and we 
infer a transition luminosity of 7 x lO^*^ erg s"', which is in- 
line w ith typical transition luminosities (Tanaka & Shibazak^ 
1 1996*). The observation 3 luminosity, which is a factor of 470 
lower than the transition flux, is only about an order of magni- 
tude above the quiescent luminosity of the black hole system 
that is brightest in quiescence (V404 Cyg, Garcia et al. 2001). 
Thus, our XTE J 1650-500 observations cover a significant 
portion of the luminosity range between the soft-to-hard state 
transition and quiescence, but we note that the quiescent lumi- 
nosity for XTE J 1650-500 could be considerably lower than 
for V404 Cyg. 

The trend observed for XTE J 1650-500 that the band- 
limited noise gradually shifts to lower frequency immediately 
after the transition to the hard state (Kalemci et al. 2003b) 
has also been observed f or several other black hole systems 
dKalemci et al 2003a: Nowak. Wilms & Dove ,2002). The 
drop in frequency can be explained by accretion models con- 



sisting of an inner optically thin corona and an outer op- 
tically thick disk. As the mass accretion rate drops, the 
inner radius of the optically thick disk (/?,„) moves away 
from the black hole, causing a drop in the dynamical time 
scales and, t hus, the characteristic frequencies for the system 
jdi Matteo&Psaltis 1999 ). After the XTE J1650-500 transi- 
tion to the hard state, the Lorentzian half-width is correlated 
with the value of the power-law index down to at least a half- 
width of 1 Hz and possibly down to 0.25 Hz (see §3.3). Cor- 
relations between frequency (either half-width or QPO fre- 
quency) and r have been seen before during RXTE observa- 
tions of Cyg X-1, GX 339-4, and other black hole systems 
in their hard states ("Gilfanov. Churazov & Revnivtsevlll 



Hevnivtsev, Gilfanov & Churazov 2001; Kalemci 2002|). As 
these authors argue, this correlation is expected for the disk 
and inner corona picture because the soft photons from the 
optically thick disk cool the corona, leading to a softer power- 
law index. Thus, as /?„, increases, and the characteristic fre- 
quencies drop, the corona is subject to a lower flux of soft pho- 
tons from the disk, and the coronal spectrum hardens. This ef- 
fectis predicted in both ADAF (Esin, McClintock & Naravaiil 
fT997l) models and in general disk and inner corona models 
such as the sphere-i-disk model (iNowak. Wilms & Dovel200^ 
where the X-ray emission from the corona is predominantly 
due to thermal Comptonization. 

While the half-width and F are correlated above 0.25-1 Hz, 
F does not get any harder than 1 .65- 1 .70 for XTE J 1650-500, 
and the spectrum softens to 1.93 ±0.13 as the half- width and 
flux continue to drop. This behavior can be interpreted within 
the ADAF or sphere-ndisk (henceforth ADAF/SD) models. 
The decrease in the half-width frequency indicates that the 
disk inner radius continues to recede; however, the disk flux 
is at such a low level that it no longer significantly cools the 
corona, and the spectrum does not continue to get harder. In- 
stead, the hardness of the spectrum depends on the "Compton 
y-parameter," which is a function of the optical depth and the 
temperature of the corona (iRvbicki & L ightman 1979). As 
the mass accretion rate drops, the coronal optical depth de- 
creases, and a softer spectrum is predicted. Softer spectra 
have been observed at low luminosities for several black hole 
sources besides XTE J 1650-500. XTE Jl 1 18H-480 shows sig- 
nificant softening that starts at an X-ray luminosity of 2 x 10 
erg s"' (ii/1.8 kpc)^, and, at the end of its 2001 outburst, 
XTE J 1550-564 begins gra dual softening a t a luminosity of 
2 X 10^^^ erg s"' (c//5.3 kpc)^ (' Kalemcil20^ . Other examples 
of spectral softening at low luminosities include GS 1 124-68 
(Ebisawaet al. 1994), as well as some b lack hole sou rces in 
or close to quiescence ( Mg^linjggk^et al. 2003; Corbe retal] 
[2003; Kong et al. 20021 iTomsick. Co rbel & Kaaret 2001^ 

Although spectral fits with the ADAF model imply rel- 
atively large changes in /?„, between the soft and hard 
states, some previous black hole timing studies conclude that 
such c hanges, which can be as large as factors of 30-3000 
llEsin. M cClintock & Naravan 1997 1), do not occur. For sev- 
eral black hole systems, di Matteo & PsaltisI ill 999 1) find that 
the characteristic frequencies do not vary by more than a fac- 
tor of 50 between states. For Keplerian orbits, /?,„ oc v~^l^, 
where v is the characteristic frequency, so that a factor of 50 
in frequency implies a change in Rj„ by a factor of 14, which 
is significantly less than indicated by ADAF fits. However, 
for XTE J 1650-500, we find that the Lorentzian half-width 
changes by a factor of nearly 1200 while the source remains 
in the hard state, implying that /?,„ can change by as much as a 
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factor of 1 13. Results of fast optical photometry for quiescent 
black hole systems suggest th at Rm may change by an even 
larger factor. For A 0620-00, Hvnes et alJ j2003l) find that the 
optical power spectrum is dominated by band limited noise 
with a break frequency of 9.5 x 10"^ Hz, which is a factor of 
3.7 lower than the lowest frequency we find for XTE J 1650- 
500. If the frequency is, in fact, relat ed to Ri„, then our re sult 
for XTE J 1650-500 and the results of lHvnes et af] ( l200 3>) im- 
ply that very large values of occur at low luminosities, 
and that the Rj„ evolution occurs gradually and continues to 
change all the way down to quiescent luminosities. These re- 
sults also imply that an extremely large value of Rj„ is not 
required for a source to show hard state properties. 

While the ADAF/SD models described above can explain 
many of the X-ray properties we see for XTE J 1650-500, 
radio observations of XTE J 1650-500 made after the tran- 
sition to the hard state indicate the likely presence of a com- 
pact jet (Corbel et al., in prep.), which is not accounted for 
in these models. Howev er, the "magnetic corona" model of 
iMerloni & FabianI (|2002) also explains several of the XTE 
J 1650-500 X-ray properties while predicting that energy 
stored in the corona should be available to drive strong out- 
flows. In this model, the X-ray emission mechanism is also 
Comptonization, so that a high flu x of disk photons leads to 
a softer energy spectrum. Thus, Mer loni & FabianI J2002h 
predict that the spectrum will harden as the mass accre- 
tion rate drops, and they also predict that the spectrum will 
soften at very low mass accretion rates as more of the coro- 
nal energy is carried away by the outflow. Recent work by 
iFender. Gal lo & Jonker' (2003) suggests that we are observ- 
ing XTE J 1650-500 in a luminosity regime where this could 
occur While t he XTE J 1 650-500 spectral evolution is con- 
sistent with the IMerloni & Fab ian (2002) model, and we also 
previously argued that this model may be able to account for 
the X-ray flares reported in 'Toms ick et alJ ( 12003 ). it is un- 
clear if the magnetic corona model can account for the fact 
that the observed characteristic frequency for XTE J 1650- 
500 changes by a factor of 1200 in the hard state. Perhaps 
this could indicate large variations in the size of the magnetic 
corona or in the sizes of the active coronal regions, but since 
Ri„ is not expected to change by a large factor in this model, 
we do not see an obvious explanation for the large frequency 
range. 

A jet-dominated model has also been suggested for black 
holes in their hard state. In this model, the compact jet con- 
tributes a large fraction of the emission from the radio band 
up to X-ray ene rgies, and a syn chrotron emission mechanism 
is invoked (Markoff. Falcke & Fender 2001) . Evidence in fa- 
vor of this model includes observations of cor related radio 
and X-ray flux for several black hole systems (' Corbel et al.l 
12003; Gallo. Fender & Poolev 2003), and good agreement be- 
tween the predictions of the model and the broadband GX 
339^ energy spectra over a large range of mass accretion 
rates iMarkoff et al. 2003). In the model, X -rays are produced 
when electrons are accelerated at a shock that is Zsh ^ ^^^Rg 
from the black hole, where Rg = GMjc^ and M is the black 
hole mass. The X-ray spectrum is predicted to have a power- 
law shape out to a cutoff energy, and the power-law index is 
set by the index of the electron energy distribution, p, where 
dNjdE cx E'P. IMarkoff etaf] (120031) find that they can ex- 
plain the radio/X-ray flux correlation for GX 339-4 by vary- 
ing only the jet power, keeping p and z,/, constant. The fact 
that r changes for XTE J1650-500 and also for GX 339-4 



llCorbel et al.l|20()l . indicates that if this model is correct, n 
likely does change with source luminosity, but we suspect that 
changes in F are probably not a major problem for this model 
as one might expect p to vary somewhat with jet power. On 
the other hand, one would expect that the characteristic fre- 
quencies of the system would be related to the size of the 
source. Assuming that the size of the X-ray emitting region is 
related to the location of the shock, Zs/,, it may be difficult for 
the jet-dominated model to explain the large range of charac- 
teristic frequencies we observe for XTE J 1650-500 without 
varying Zjft. 

5. CONCLUSIONS 

We have compared the emission properties of XTE J 1650- 
500 in the hard state over a large range of luminosities to 
ADAF/SD, magnetic corona, and jet-dominated models. Al- 
though the ADAF/SD model provides natural explanations for 
the X-ray spectral and timing properties, these models do not 
include outflows and thus cannot provide a complete descrip- 
tion of the system. In addition, ADAFs may be convectively 
unstable as mentioned above. The magnetic corona model, 
which does incorporate a jet, predicts a similar spectral evo- 
lution to the ADAF/SD model, but it it unclear if it can ex- 
plain the large change in the Lorentzian half-width that we 
see for XTE J1650-500. The jet-dominated model can ex- 
plain X-ray/radio correlations that are seen for several black 
hole sources without varying Zs/,; however, it is doubtful that 
the model can reproduce the observed X-ray timing properties 
with a constant value of z,/, . 

Regardless of which model is correct, the fact that the gen- 
eral shape of the energy and power spectra remain the same 
while the parameters evolve indicates that the X-ray emitting 
region of the system likely keeps the same overall structure 
even though the parameter changes indicate that some physi- 
cal properties change, such as the inner radius of the accretion 
disk, the size of the corona, or the shock location, consider- 
ing the models described above. It is clear that the source 
retains its hard state X-ray properties down to a luminosity 
of 1.5 X 10^"^ erg s"', and we do not see any evidence for a 
sharp transition to the quiescent state. Furthermore, our re- 
sults for XTE J 1650-500 combined with the determination 
that the optical power spectrum for A 0620-00 in quiescence 
exhibits band li mited noise w ith an even lower characteristic 
frequency (Hv nes et alJl2003l) suggests that black hole tran- 
sients keep the same overall structure into quiescence. Thus, 
these results favor the hypothesis that the quiescent state is 
simply a low luminosity version of the hard state. 
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APPENDIX A 

We performed simulations to study the effects of pile-up on the observation 3 ACIS power spectrum. The simulations have 3 
main goals: I. To understand why the Poisson noise level for the actual Leahy normalized ACIS power spectrum is 1.77; 2. To 
determine how much the fractional rms amplitude is reduced by pile-up; and 3. To check that pile-up does not alter the overall 
shape of the power spectrum. Our goals were accomplished via 3 simulations. 

Simulation #1: We produced a simulated Poisson noise light curve with the same time resolution as the actual data (but with 
much longer duration to improve the statistics). According to the binomial probabilities, we produced a piled-up light curve. The 
2 adjustable parameters are the input count rate (before pile-up) and the probability that any pair of events occurring in the same 
time bin will pile-up. We adjusted these parameters to match the actual ACIS observation as shown in Figure0 We produced a 
Leahy normalized power spectrum using the simulated light curve, and the measured Poisson level is L78. 

Simulation #2: In addition to the deadtime caused by pile-up, there is a fixed, instrumental deadtime that is independent of the 
event rate and is set by the CCD readout time of 4L04 milliseconds. To determine the effect of the readout time, we simulated a 
light curve where there is a probability (equal to the readout time divided by the 0.4 s exposure time per frame) that a photon that 
hits the detector will not be counted. We produced a Leahy normalized power spectrum using the simulated light curve, and the 
Poisson level is L97. 

Simulation #3: Here, we made simulated light curves similar to those used for simulation #1 except that we included white 
noise in excess of the Poisson noise. We produced a light curve free of pile-up and also a light curve with pile-up and made rms 
normalized power spectra for each. The fractional rms amplitude changes from 19% to 18% when pile-up effects are included. 
The shapes of the power spectra with and without pile-up are not significantly different. 

In summary, these simulations demonstrate that the low Poisson noise level measured for the actual data is mainly due to 
pile-up effects, and there is also a small drop due to the CCD readout time. For noise levels similar to those we measure for 
observation 3 (~18%), pile-up only causes the fractional rms amplitude to drop by about 1%, and the pile-up does not change 
the shape of the power spectrum. 

APPENDIX B 

Here, we describe how we corrected the observation 3 ACIS response matrix for the energy dependence of the Chandra PSF. 
We generated this correction using the web-based Chandra calibration tool "ChaRT" (Chandra Ray Tracing). For 13 different 
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Fig. 7. — The solid histogram shows the fraction of 0.4 s observation 3 light curve bins with with 0, 1,2, or 3 counts/bin. The dashed histogram shows the 
same for a simulated light curve that includes the effects of photon pile-up. This illustrates that our simulated light curves are a good match to the actual data. 
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Fig. 8. — ChaRT (Chandra Ray Tracing) results for a source 2'. 7 off-axis showing the fraction of counts in the elliptical annulus we used to extract the 
observation 3 ACIS spectrum (see text for the precise region parameters) as a function of energy. The points come from directly from the ChaRT simulations, 
and the solid line is the spline interpolation that was applied to the ACIS response matrix for observation 3 to account for the energy dependence of the Chandra 
PSR 



energies from 0.2 to 10 keV, we used ChaRT to simulate a monochromatic source at the position of XTE J1650-500, 2'. 7 off- 
axis. For each energy, we determined the fraction of counts in the same elliptical annulus that we used to extract the XTE 
J 1650-500 spectrum (see text), and the fractions are shown in Figure |8l We then used a spline interpolation to determine the 
fraction for each energy bin in the response matrix, and we applied the interpolated values to the response matrix. While this 
correction is necessary, we note that, in this case, it only changes the measurement of the power-law index from r = 1.79±0.13 
to r = 1.93 ±0.13. Although we do not have a method of accurately determining the systematic error on this correction, even if 
the correction is only good to 10-20%, the systematic error is much smaller than the statistical error 
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Table 1 

Chandra and KXTE Observations 



Obs./ 






Energy Bwd 


Exposure 


Count 


UTDate 


MID Start" 


Instrument 


(keV) 


(s) 


Rate* 


1 


52298.00688 


ACIS/HETG/MEG 


0.7-7 


10001 


0.479 ±0.007 


(2002 Jan. 23-24) 


52298.00688 


ACIS/HETG/HEG 


1-9 


10001 


0.271 ±0.005 




52298.03298 


RXTEI?Ch 


3-20 


1904 


9.42 ±0.14 


2 


52309.62048 


ACIS/HETG/MEG 


0.7-7 


9509 


0.457 ± 0.007 


(2002 Feb. 4) 


52309.62048 


ACIS/HETG/HEG 


1-9 


9509 


0.248 ±0.005 




52309.71204 


RXTE/PCA 


3-20 


1872 


8.69 ±0.14 


3 


52335.10444 


ACTS'' 


0.5-8 


1 8266 


0.493 ±0,005'' 


(2()()2 Mar. 2) 


52335.15312 


KXJE/PCA 


3-20 


2384 


1.38 ±0.11 



"Modified Julian Date (JD-240()0()0.5) at exposure midpoint. 

''The MEG and HEG rates quoted are for the -l-l and —1 grating orders combined. 

■^For this observation, we used a 1/8 subarray and an offset pointing of 2'. 7 to mitigate the effects of photon pile-up. 

''This is the count rate for the source, including the core, which is moderately affected by photon pile-up. We estimate that the count rate would be ~15% higher 
without pile-up. 



Table 2 

Parameters from Chandra and PCA Spectral Fits 





Nh 




Flux 




Fe Ka Line EW 


Obs." 


(10^' cm--) 


r 


(1-9 IceV)* 




Upper Limit (eV)*^ 


1 


6.7 ±0.5 


1.66 ±0.05 


4.76 X Ur" 


116/155 


< 162 


2 


6.7 ±0.5 


1.66 ±0.05 


4.46 X ur" 


164/155 


< 169 


3 


6.4 ±0.7 


1.93±0.13 


7.9 X 


104/100 


<212 


3 


6.7'' 


1.96 ±0.09 


8.0 X 10-'2 


104/101 


<212 



"These are parameters from absorbed power-law fits to the spectra. AH errors and upper limits are 90% confidence (Ax^ = 2.7). 
''The 1-9 keV absorbed flux in units of erg cm~^ s~' . 

^^90% confidence upper limit on the equivalent width of a narrow emission line at 6.4 keV. 
''Fixed. 



Table 3 

Paramaters for the Power Spectra 



Obs." 


Instrument 


Energy Band (keV) 


FWHM/2'' 


Fractional RMS (%) 




1-h2'^ 


PCA 


3-20 


0.067 ±0.007 


33.1 ±1.1 


93/72 


1-H2 


MEG-hHEG 


0.7-9 


0.079 ±0.020 


30.1 ±2.5 


25/23 


3 


acis 


0.5-8 


0.0035 ±0.0010 


19+^ 


56/65 



"These are parameters from fits with a zero-centered Lorentzian. Errors are 68% confidence. 
''The half-width of the zero-centered Lorentzian in Hz. 

■^A QPO may be present near 0.1 1 Hz with a fractional rms of about 12%. The significance of the QPO is only required at a confidence level of between 2-a 
and 3-a. 



